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We propose a general description on the unambiguous discrimination of mixed states according to 
the system-environment couphng, and present a procedure to reduce this to a standard semidefinite 
programming problem. In the two states case, we introduce the canonical vectors and partly simplify 
the problem to the case of discrimination between pairs of canonical vectors. By considering the 
positivity of the two by two matrices, we obtain a series of new upper bounds of the total success 
probability which depends on both the prior probabilities and specific state structures. 
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I. INTRODUCTION 

Quantum state discrimination (QSD) is one of the 
fundamentally important problems in quantum informa- 
tion science. Especially in quantum communication and 
quantum cryptography, many novel schemes are based 
on the fact that nonorthogonal states cannot be discrim- 
inated determinately. Henceforth study on the discrimi- 
nation of quantum states has a close relation to the secu- 
rity of quantum cryptographic protocols. On the other 
hand, since there is no measurement that can perform 
a perfect identification, several strategies have been pro- 
posed in QSD based on different criteria. One of these 
is the minimum-error discrimination which permits 
incorrect outcomes during the measurement procedure. 
The other one is unambiguous discrimination (UD) of 
quantum states. This sort of discrimination procedure 
never gives an erroneous result, but sometimes it fails. 
Here we consider the latter case which has received much 
attention recently. 

In thepure states case, UD has been widely consid- 
ered 0, Q • While in mixed states case, it seems to be 
a hard problem. In many earlier works, some useful 
bounds of the total success probability P, together with 
several useful reduction theorem, have been presented 
i, S, i, 0, i, i, El- However, totally solving this 
problem seems not so easy at all. What's more, even 
for the simple case, e.g., UD between two mixed states, 
which has been wildly studied recently, there still exists 
many questions which are not clear to us. 

The standard description of UD among mixed states 
are usually formulated as this: given a set of mixed 
states {pi, P2, ■ ■ • , PAf} with the corresponding prior 
probabilities {771, 772, . . . , ?7Ar}, the aim of discriminat- 
ing these states unambiguously is to find a iV + 1- 
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element positive operator value measurement (POVM) 
{So, E^,E2, ...,En} with E, = I and TT{EkPi) = 

Pk5ki{k,l ^ 0), such that the measurement operator Ek 
gets a result with the probability pk only when the in- 
put state is pk- Here Eq denotes the inconclusive mea- 
surement where the identification fails. The average 
failure probability is described hy Q = X^s^i Qk with 
Qk = ?7fcTr(i?o/Ofc) being the failure probability of identi- 
fying Pk. Equivalently, one can also concentrate on the 
total success probability P = 1 — Q = J2k=i VkT^T^{EkPk)- 

From the general viewpoint, UD can be regarded as 
some kind of physically accessible transformation on a 
finite number of input states ^ : {pi, p2, ■ ■ ■ , Pn} 
{(Ti, (72, ■ • ■ , Cat} The major character of this trans- 
formation is that it is probabilistic, accurate, and the 
output states ak must be orthogonal to each other so 
that they can be identified perfectly. There are several 
equivalent approaches to describe a completely positive 
(CP) map [13]. For example, it can be represented in a 
Kraus operator sum form. Also it can be implemented by 
employing a unitary transformation on the system plus 
ancilla, i.e. ^(p) — TrE'[Up(^ peW^ I ® Pe'], where pE 
is the initial state of the ancilla system, / denotes the 
identity operator in output Hilbert space 7^2, Pe' is a 
projector \tiT-Le', and TCii^He ='H2®'He'- 

In this paper, we consider to discriminate mixed states 
unambiguously from the system-ancilla coupling view- 
point. By constructing the whole unitary transformation 
on the combinations of the inputs and the auxiliary sys- 
tem, we obtain the necessary and sufficient conditions on 
the existence of a UD strategy. We point out that in the 
more general case to find the optimal UD strategy can be 
reduced to a standard semidefinite programming prob- 
lem. Especially, in the case of UD between two mixed 
states, we obtain a series of new upper bounds of the 
success probability which are closely related to the struc- 
ture of the input quantum states together with the ratio 
of prior probabilities. In some sense our result confirms 
the conjecture made by Bergou et al [Tlj . 
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II. GENERAL DESCRIPTION OF 
UNAMBIGUOUS DISCRIMINATION 

Let us start with the general UD of N mixed states. 
For any mixed state pk, it can always be regarded as the 
mixture of pure states, i.e. pk = J2m IV'm^) (V'™'' L where 
I'i/^m'') are nonnormalized state vectors. Here for sim- 
plicity, we assume that \'4>m^) are linearly independent. 
Firstly, we suppose the intersection of the supports of 
two density matrices pk and pi is empty (except for a 
trivial zero vector). This also indicates that all vectors 

{IV'm''): ■ • ■ I IV'n'')} are linearly independent. By intro- 
ducing suitable auxiliary system, we consider the follow- 
ing unitary realization of the CP map ^ 



f/|^«)l|O)^ = |0W)2a|Po)p + |/3L')) 



2ap- 



(1) 



Here He' = 'Ha ® "Hp, |0) is the fixed initial state of the 
environment, |Po) is the state of the probe system satis- 



tying (Pol A 



0, and we also use the tilde to denote 



a nonnormalized state vector. The output state can 
be obtained by tracing over the subsystem a after we get 
a measurement outcome corresponding to the probe \Po), 
i-e- afc=E™Tr<,[|<>L'^)(0L''|]. 

On the other hand, if the intersection of the supports 
of two density matrix pk and pi is not empty, there exists 
at least one state vector g supp(pk)C\ supp{pi). From 
the definition of the CP map, we have 



C/|V)|0) 



|<^('-))|Po> + |/3('=') 
|<^('^)|Po) + |/3(')). 



(2) 



Since 10'^'"^) ^ 10*^'^) (The output states are differ- 
ent from each other), Eq. ^ is satisfied only when 
\(p^^^) — l^^'-*) = 0, hence any state contained in 
supp{pk) n supp{pi) has no contribution to the desired 
transformation. Thus it's enough to consider the case 
supp{pk)C\ supp{pi) = {0}, which reproduces the known 
results Q. 

The inner-product preservation of unitary transforma- 
tion leads us to the following equation 



X = B > 



(3) 



with 



w 



I Wkk ... Wkl 



\ wik ■■■ Wll 



{we{X,Y,B)}. (4) 



Here liki are all block matrices with {Xki)mn — 
(V^^^l^i'^ {YklUn = i&^n^), and (BklUn = 
(/3m "* |/3n "*) respectively. Also we can find that all the three 
matrices (X, Y , B) are Hermitian, and positive semidef- 
inite. Since Ufc are orthogonal to each other, we have 
{'pm}\4>7i') = 0(fc I). This indicates Y is quasi-diagonal 
and can be written as 1" = diag{Ykk, ■ ■ ■ ,Yii}. 



Contrarily, if there exists a positive semidefinite Y ma- 
trix satisfying Eq. ([3]), we can always choose suitable 

state vectors \4>rn) and |/3m'') such that X — B + Y. 
With the standard Gram-Schmidt orthogonalization pro- 
cedure, the desired the unitary transformation can be 
easily obtained. 

We conclude the above discussion by the following the- 
orem. 

Theorem 1. N mixed states {pi, p2, ■ ■ ■ Pn} can be un- 
ambiguously discriminated if and only if there exists a 
positive semidefinite quasi- diagonal matrix Y such that 
X ~ Y > 0. Moreover, if the input states are chosen 
with prior probabilities {?7i, ?72, • • • w} i^*^^ l^kVk — 1; 
the total success probability will be P = J2k Vk Tr(Ykk)- 

This theorem characterize the general properties of UD 
among N mixed states in the system-ancilla framework. 
One can also easily check that it is consistent with ear- 
lier works 0, H, [a 0, IE i, 0, [HI- In a more realistic 
situation, people often concentrate on the total success 
probability of such physical transformation. This indi- 
cates that we should make the probability P as high as 
possible. Mathematically, this is equivalent to maximiz- 
ing 



P : 



y^TyfcTr(yfcfc) 

k 



(5) 



under the constraints 



X -Y >0,andY >0. 



(6) 

Usually given the input mixed states, we can get to know 
the matrix X exactly. Therefore the only thing we should 
do is to find the optimal positive semidefinite matrix Y 
which maximizes the success probability P. By redefin- 
ing a series of new matrices Fq = diag{X,0}, F^*^ — 
diag{EP'',-EP'^}, and Gf = dmff{zFf , -iFf }, where 
i is the basic imaginary unit, and -E^* are matrices corre- 
sponding the block matrices Ykk with {E^'^)mn = SmpSnq, 
the problem under consideration can be reformulated as 

^VkMYkk) (7) 



max 
Y 



k 



subject to 

^^0 



E 

k.p,q 



{Remk)pq]Fj:'^ + Im[{Ykk)pq]Gf} 



>0, 



where Re[{Ykk)pq] and Im[(Ykk)pq] represent the real 
and imaginary parts of the matrices elements {Ykk)pq 
respectively. This is a standard semi-definite program- 
ming (SDP) problem ^14] , and can be solved by numeric 
method efficiently (one can also find another method to 
reduce this to a SDP problem in TF| , which is equivalent 
to our result). Therefore in principle, the optimal success 
probability of UD of mixed states can be found numeri- 
cally. Actually once we have found the optimal matrix Y , 
with the standard procedure, we can construct the cor- 
responding unitary implementation of the discrimination 
operation. 
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III. UNAMBIGUOUS DISCRIMINATION OF 
TWO MIXED STATES 

In the above discussion, we have given a general de- 
scription on UD among N mixed input states. To be 
specific, in the following, we will focus on a particular 
case, i.e. UD between two mixed states. This is a ba- 
sic and very important case in the study of UD, and 
much attention has been paid to the problem recently. 
In Rudolph et al. present the lower bound on the 
failure probability Q, and later, it has been pointed out 
that there exist mixed states for which the lower bound 
can not be reached for any prio r probabilities. Based on 
these facts, Raynal et al. investigated a large class 

of two mixed states discrimination, they also found the 
necessary and sufficient conditions for two mixed state 
to saturate these bounds. In all these works, Q is con- 
sidered in three different regions, depending on the ra- 
tio between two prior probabilities. Recently, Bergou et 
al. [ll| have considered the discrimination of two sub- 
spaces and they find that for this special case there are 
many parameter regions which can give different mini- 
mal failure probabilities. The regions depend on both 
the prior probabilities and the specific structure of the 
two subspaces. The lower bound 2 ^771 772 F of the failure 
probability Q can be reached only when the prior proba- 
bilities lie in some specific regions. Later they conjecture 
that this phenomenon occurs for any two mixed states. 
In the following, we will show that this result is indeed 
universal. 

When restricted to two-state case, Eq. ([3]) can be sim- 
plified as 

where X^,; arise from the decompositions of pi 
and p2- Usually there exist many other en- 
sembles which can generate the same operators, 
i.e. PI = = 

(ii/;j^)),iv;«),...)[/tc/((^(i)|,(^(i)|,...r 

(|fi),|r2),...)((fi|,(f2|,...f , where U {V for p^) 
is a unitary matrix and T represents the transpose of 
the matrix. This is known as the unitary freedom for 
density matrices [l^. Hence we can also write down 
the correspondence of Eq. ([3]) according to this new 
decomposition 

X'-Y'>Q. (9) 

Since X' = diag{?7, T^}Xdiag{[/^, y^}, we can immedi- 
ately obtain that this will not affect the total success 
probability P we consider here (This is also general for 
N input mixed states). 

Keeping in mind that U and V can be arbitrary, we 
can choose the two matrices appropriately such that 
UXi2V^ = diag{diag{/i,/2,...,/t}, 0}, where we as- 
sume Xii and X22 are u x u and v x v matrices respec- 
tively with u < V, fm are the singular values of X12, and 



is a {u — t) X {v — t) zero matrix. This implies there 
exist some kinds of decompositions of pi and p2, namely, 
Pi = J2m \rrn){fm\ and p2 = |sn)(s„|, which satisfies 
the following equations 

/f Is \ — J fm^rnn {m,n) < t, , . 

Krm\Sn) - I Q otherwise. 

The singular values fm have very interesting properties 
and we characterize this by the following theorem [17]. 

Theorem 2. Given two mixed states density matri- 
ces pi and p2, there exist two sets of canonical vectors 
{1^1)7 1^2), ■ • ■} OJ^c? {|si), \s2), ■ ■ ■}, which generate pi and 
P2 respectively, such that Eq. \10\} is satisfied. And the 
fidelity of the two density matrices can be formulated as 

F = Tr^J p\''^P2p\'^ = Em /m- 

Proof: The only thing we should do now is to prove 
the second part of this theorem. Consider the spectral 
decompositions of pi and p2 

pi=Y,a,\a,){a,l P2 = ^ (/3/.|. (11) 

i h 

According to the definition of fidelity _F, we obtain 

■^^{f3h\aj). (12) 

Now we definite a new matrix 

A^h = Va'^^/fh{a^\(3h). (13) 
The fidelity F can be rewritten as this 

F = Tr^JpY^p2pY^ = TrVHt. (14) 

Since ^ is a complex matrix, using a sin- 
gular value decomposition, we have A ~ 
[/idiag{diag{/i,/2,...,/t},"0}Fi with fm > for 
all 1 < m < Thus the fidelity becomes 

F = TrVZ4t =^/,„. (15) 

m 

On the other hand, because of the unitary freedom in the 
ensemble representation of density matrices, we can find 
two unitary operations U2 and V2 such that C/2-'^i2V2 = 
A. Therefore A and X12 have the same singular values, 
which completes the proof. 

Theorem 2 indicates for any two mixed states, it is 
always possible to find two sets of canonical vectors 
{|fi), |f2), . . .} and {|si), |s2), . . .} so that |f,„) only have 
a nonzero overlap with |s„i)- When (n,m) > t, one can 
easily check that \ fm) and |s„) lie in the subspace orthog- 
onal to the supports ofp2 and pi respectively. From the 
reduction theorem in [7|, we conclude that UD between 
pi and p2 is equivalent to that between the two newly 
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defined density matrices p\ = X]L=i 1^™) ^^"^ 

P'l = ELi |Sn)(Sn|/iV2 with Nx = Tr(X;L=l l^mX^ml) 

and Ni = Tr(^yj^-^ |s„)(5n|) being the corresponding 
normahzation factors. According to the system-ancilla 
model (Theorem 1), Eq. ([5]) can always be reduced to a 
2t X 2t matrix 



Xn-Fn diag{/i,...,/t} 
diag{/i, . . . 



> 0, (16) 



where we have used the same notations for simplicity. 

Equation (|16p supplies enough information which can 
be used to demonstrate our main results. Actually, since 
X — y is positive semidefinite, from the standard linear 
algebra theory, we have that every principal minor of 
X — Y \s also positive semidefinite, i.e. 



fn 



> 0, 



(17) 



where we have used y„i and Zm to denote the diagonal 
elements of Yn and I22 respectively, and (sm) is the 
modulus of the vector |fm) (|sm))- Therefore by intro- 
ducing the canonical state vectors, the question can be in 
part reduced to UD between pairs of state vectors \ fm) 
and \srn)- Such question has been solved in many earlier 
works, and the results are listed as follows 

^ Tjmax 
— 



(18) 



The above expression shows that for every m, the 
maximal value that Pm can achieve has a close rela- 
tion with the specific configuration of y^^, i"m, Sm, and 

fm- Generally, for different m, P,„ will have very differ- 
ent expressions. Therefore, the total success probability 
P = J2m ^™ cannot always be represented as a function 
of the fidelity F = S„/™. 

Specifically, in the following we will focus on some 
special cases. Firstly, if for all m = l,...,t, we have 



^2—7^1 then according to the above equation, 

the upper bound of the total success probability can be 
rewritten as 



p = E^™^E^^ 



= 2^ Virm + ?72S„i - 2 ^771 772 /m 

rn 

= 1-20^F. (19) 

The corresponding lower bound of the failure probability 
becomes Q = 1 — P > 2 .^771 772 -F. This bound has been 



proved to be the minimal value of Q for any type of in- 
put configurations. However, our result shows that even 
in this special case, the lower bound of Q can only be 
possibly saturated. This occurs, for example, when the 
canonical vectors are orthogonal to each other. In gen- 
eral case, since Xn and X22 are not diagonal matrices, 
this lower bound cannot always be reached. 

In the second example, we assume that > 

for all 1 < m < A simple algebra will lead us 
to the following bound of the total success probability 
P ^ 7?i(l — 12m fm/ ^m)- If wc introducc a new op- 
erator C2 = J2m \sm){sm.\ composcd of thc normalized 
canonical vectors of p2, we can reformulate P as P < 
7yi(l - Tr(piC2)), or equivalently Q > 772 + 77iTr(piC2). 
When \sm) are orthogonal to each other, C2 is nothing 
but the projection onto the support of p2- 

Thirdly, if we have < 7^ (V777,) , the total success 

probability satisfies P < 7^2(1 - = '72(1 - 

Tr(p2C'i)) with Ci = J2„i\'^rn){rm\- Correspondingly, 
the failure probability becomes Q > rji + 772Tr(p2Ci)- 

For mixed states pi and p2, we always have (r^, Sm) < 
1. This indicates that the failure probability Q can 
never reach the bound 2 ^771 772 F for the latter two cases. 
Generally, different canonical vectors of the input states 
will separate the parameter space into different regions, 
and the lower bound of Q is determined by both the 
prior probabilities and the structure of states. More- 
over, in each region, the lower bound of Q can not al- 
ways be reached. Mathematically, to judge whether the 
lower bound can be saturated is equivalent to determin- 
ing whether there exists a positive semidefinite matrix 
F > such that Eq. ([8]) is satisfied. This problem is 
often called semidefinite feasibility problem (SDFP). Un- 
fortunately, the complexity of SDFP is still not known, 
and currently we can only say that it cannot be a NP- 
complete problem unless NP=NP-complete. Therefore 
to judge whether the bound of Q can be reached or not 
seems to be a hard problem. But in some special case (for 
example, the canonical vectors are orthogonal to each 
other, or AT — y is a diagonally dominant matrix), some 
known results in linear algebra theory will be helpful to 
solve this problem. 



IV. EXAMPLES 

In many related works, the upper bound of thc success 
probability P is only considered in three different inter- 
vals, which depends on the ratio of 771 and 7^2 together 
with the fidelity F and supports of the input states. Here 
by introducing the decomposition in Theorem 2, we find 
a series of parameter regions related to the specific in- 
put states. In addition, from the system-ancilla coupling 
viewpoint, one can also derive the corresponding results 
in [95 and For example, if we definite B = X' — Y' 
in Eq. then since B is positive semidefinite, we have 



5 



\/»7l/^2 


P 


PRa{PRu) 


[0, cosei] 


i772(sm^6li +sm^6'2) 


ir72(sm^6li +sm^6l2) + 

2^1 COS-^fl'+COS-4, (O'' ^ )) 


[cosei,Rai{F)] 


^ - ^JVin^cosOi + ir72Siri^6'2 


[Rai{F),cos92] 


i — ^^rjiri'zr ^or i — z,^r]ir)2r ) 


[cosSa, cose J 


1 - 2y/rjrn^F 






[^«2(l/i^),cok] 


i»7i(sm^6li +sin^6l2) + 
W^SI^^^^€ (or mil -F^)) 


[cosei ' 


i?7i(sin^6li +sin^6'2) 



TABLE I: Bounds of the maximal success probabilities presented in several related works. Here P denotes the bound according 
to Eq. (|18p . PRa and Pru are the results obtained from Ql and respectively. F — (cosSi +cos62)/2 is the fidelity of the two 
input mixed states. Rai = Tr(Pip2)/-F — (cos^6'i + cos^S2)/(cos0i +C0SS2) and Ra2 = -F/Tr(P2Pi) = (cos^i + cos02)/(cos^6'i + 
cos'^^2) are parameters according to [3] with Pi and P2 being the supports of pi and p2 separately. 



^Tr(i3ii)Tr(S22) > |Tr(Bi2)| for any kind of decompo- 
sitions of pi and p2 (for the definitions of Bij, see Eq. 
©). Therefore we obtain y'Tr(i?ii)Tr(B22) > where 

equality holds only when Bn — ai?22 with a € R. This 
also indicates that the output states corresponding to the 
failure measurement results cannot be used for further 
discrimination operations, which is consistent with the 
discussions in Q. To reveal the relation and difference 
between the bounds list above and those in the previ- 
ous works, in the following we will investigate a specific 
example. 

Consider two rank-2 mixed states pi — ^{\ri){ri \ + 
|r2)(r2|) and p2 = ^{\si){si\ + |s2)(s2|) with (ri|s2) = 



(?'2|si) = 0, (ri|si) = cos^i, and (r2|s2) = cos6'2. To 
simplify our consideration, we also assume (r'i|r2) = 
(si|s2) = 0. Actually, discrimination of such kind of 
mixed states has be extensively studied in Here 
we also use it to manifest the difference of the upper 
bounds presented in several related works. Suppose 
< cos^i < COS02 < 1- Then based on our former dis- 
cussions, optimal success probability P can be obtained 
exactly in five parameter regions, [0, cos^i], [cos6'i, cos6'2], 
[cos6'2, l/cos02], [l/cos02, l/cos^i], [l/cos^i , cx)) . Alter- 
natively, one can also obtain the corresponding upper 
bound of P according to [9|] and |5|]. Table 1 shows the 
details of P in each of these regions. 



The above table shows that when cos6'i — cos6'2, 
the three bounds P, Pru, and Pru are equal to each 
other. However, for the general case 008611 cos02, 
one can easily obtain P < Prq and P <__Pr«, 
and equalities hold only when cos02 < 1/771/772 < 
l/cos02- For example, if cos^i < \/'r]\lr\i < 
Rai, we have Prq — P ^ 5772 [2a;cos6'i — 008^611 — 
2x^cos0icos02/(cos^0i + 008^^2)] = f{x), where we 
have assumed x = \J 771/772. Since /(cos^i) — 
ri2COs'^9i{cos9i - cos6'2)^/[2(cos26'i + cos26'2)] > and 
f{Rai) = ?72Cos26'i(cos6'i - cos6l2)V[2(cos6li + cos6'2)^] > 
0, one immediately sees that P < P^a for any \Jr\\lr\2 m 
[cos6'i, i?ai]. These observations indicate that the bound 
presented in this work is independent of those of former 
works, and sometimes it can provide tighter bound of the 
total success probability P, as we have expected. 



and presented a procedure to reduce this to a standard 
SDP problem, which makes the problem to be solvable 
numerically. On the UD between two mixed states, we 
have introduced the canonical vectors and partly reduced 
the original problem to the UD between pairs of canoni- 
cal vectors. We present a series of new upper bounds on 
the total success probability which depends on both the 
ratio of the prior probabilities and the input state struc- 
tures. This indicates that the results in are universal 
for any type of input states. It also should be mentioned 
that throughout the paper we mainly concentrate on the 
diagonal elements of the corresponding matrices. In prac- 
tice, the non-diagonal elements also play important roles 
which deserves further investigation. 
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